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Abstract 

Classification of dark energy models in the plane of w and w' , where w is the dark energy equation 
of state and w' its time-derivative in units of the Hubble time, has been studied in the literature. 
We take the current SN la, CMB and BAO data, invoke a widely used parametrization of the 
dark energy equation of state, and obtain the constraints on the w~w' plane. We find that dark 
energy models including the cosmological constant, phantom, non-phantom barotropic fiuids, and 
monotonic up-rolling quintessence are ruled out at the 68.3% confidence level based on the current 
observational data. Down-rolling quintessence, including the thawing and the freezing models, is 
consistent with the current observations. All the above-mentioned models are still consistent with 
the data at the 95.4% confidence level. 



* Electronic address: If90222025@ntu.edu.twl 



1 



I. INTRODUCTION 



Compelling evidence from different types of observation shows that the expansion of the 
universe is accelerating at late times (see [l| for a review). Within the framework of general 
relativity, this indicates that there should exist an energy source with significant negative 
pressure, termed dark energy, to drive this acceleration. The nature of dark energy is 
generally regarded as one of the most tantalizingproblems in cosmology. Many dark energy 
models have been proposed and studied (see Q], and references therein). While the 
cosmological constant remains the simplest realization of dark energy, current observations 
do not rule out the possibility of time-evolving dark energy 

In the pursuit of revealing the nature of dark energy, cosmological observations serve to 
constrain the behavior of dark energy. Theoretical studies, on the other hand, should deter- 
mine whether dark energy models can be distinguished by their observational consequences. 
The ratio of pressure to energy density for dark energy, the equation of state w = p/p, is 
the characteristic of how the energy density evolves with time. The cosmological constant 
relates to the constant equation of state w = —1, while other dark energy models generally 
have time-evolving w. The time-derivative of w in units of the Hubble time, w' = dw/dlna, 
characterizes the dynamical behavior of the equation of state. Studies of the dynamical be- 
haviors and classification of dark energy models in the w-w' phase plane have been carried 
out J]- 8|. It is found that different dark energy models are bounded in different sectors in 
the w-w' plane. 

In this paper, on the one hand, we gather the bounds for various dark energy mod- 
els in the w-w' plane. On the other hand, we obtain the constraints on the w-w'plane 
in the redshift region < ;z < 1, by adopting a widely used parametrization [s], O, 
w{z) = wq + Wa{l — a) = wq + Wazjil + z) , bascd on the current observational data. The 
data set we use includes the recently compiled "Constitution set" of Type la supernovae 



id, 



(SN la) data [11|]-[17|, the cosmic microwave background (CMB) measurement from the 
five-year Wilkinson Microwave Anisotropy Probe (WMAP) p^, and the baryon acoustic 
oscillation (BAG) measurement from the Sloan Digital Sky Survey (SDSS) [3| and the 2dF 
Galaxy Redshift Survey (2dFGRS) j2o|. We then compare the dark energy models with the 
constraints on the w-w' plane for < 2; < 1. The work close to ours is that of Barger et 
al. 21| , in which they used the earlier data set and examined the dark energy models in the 
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WQ-Wa plane only at the redshift z = 1 . 



II. CLASSIFICATION OF DARK ENERGY MODELS 



Quintessence 



The quintessence model 22|]-2J|, which invokes a time- varying scalar field, generally 
allows its energy density and equation of state to evolve with time, and has w > —1. The 
equation of motion of the quintessence field is (p + ?>H(j) + \^<^ = 0, where H = a/a is the 
Hubble expansion rate, and V^j, = dV/dcj). In terms of w and w', the equation of motion can 



be written as 



25| 
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where the minus sign corresponds to > and the plus sign to the opposite, ^^(a) is 
the dimensionless energy density of the quintessence field, and Xq = (1 + w)/(l — w). For 
the down- rolling quintessence field (V < 0), the left-hand side of Ec|^ is positive, and 
the bound of w and w' can be obtained as w' > —3(1 — w){l + w) 



6|. The up-rolling 



quintessence field {V > 
the tracker quintessence 



3) takes the other side, w' < —3(1 — w){l + w). The bound of 
H i= obtained in flg Howeve, st.ong acceleration today, with 
w < —0.7, requires the breakdown of tracking [71]. The bound should only apply to the 
high redshift [7i], z ^ 1, which is not the region of interest in this paper. A conjectured 
limit of quintessence has been proposed in [j] as V/{—V^) < Mp, where Mp is the Plank 
mass. However, the physical origin of this limit is not clear [3] . We therefore do not impose 
this constraint on the quintessence model. Caldwell and hinder identified two categories of 

n 

quintessence models, "thawing" and "freezing" , based on their dynamical behavior [4] . For 
the thawing models, the equation of state is ~ — 1 at early times, but grows less negative 
with time as w' > 0. The bounds of the thawing models are {1 + w) < w' < 3(1 + w). For 
the freezing models, initially the equation of state is w > —1 with w' < 0, but the field is 
frozen at late times where w ^ —1 and w' ^ 0. The bounds of the freezing models are 
3w{l + w) < w' < 0.2w{l + w). Note that the upper bound for the freezing models is only 
valid for z < 1. 
Phantom 

The phantom model has negative kinetic energy and the equation of state w < —1 
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The equation of motion of the phantom field is + 3iJ0 — = 0. In terms of w and w', 
the equation of motion can be written as 



± 
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-3(1 +«;) 
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1 (iln(xp) 
6 dln{a) 



(2) 



where the plus sign corresponds to > and the minus sign to the opposite, f^^(a) is the 
dimensionless energy density of the phantom field, and Xp = —(1 + w)/{l — w). For the 
up- rolling phantom field (V > 0), the left-hand side of Eq. ([2]) is positive, and the bound 
of w and w' can be obtained as w' < —3(1 — 'w){l + w). The down-rolling phantom field 
{V < 0) takes the other side w' > —3(1 — w){l + w). Note that Eq. ([2]) and the bounds are 
different from those obtained in Q| . 
Barotropic fluids 

Barotropic fluids are those for which the pressure is an explicit function of the energy 
density, p = f{p) (see |8|] and references therein). The expression for w' can be written 
as [s], Ql 



w 



-3(1 +w) 



dp 
dp 



w 



(3) 



The sound speed for a barotropic fluid is given by = dp/ dp. To ensure stability, we 
must have > 0, which gives the bound w' < 3w{l + w) for non-phantom {w > — 1) 
barotropic fluids |5|, l8|]. For causality, we further require cj. < 1 29|, which gives the bound 
w' > -3(1 +w){l- w) for > -1 y. 

The classification of the above-mentioned dark energy models in the w-w' plane is shown 
in Fig. [H Note that all of the bounds are valid at late times for < 2; < 1. 



III. CONSTRAINTS ON THE w-w' PLANE 



A. Observational data 



We use the combined data set from three types of observations including the SN la 
observation, the CMB measurement, and the BAO measurement. We assume that the 
universe is flat in this paper. 



We use the Constitution set of SN la data compiled by Hicken et al. 



Il|-(l7|, which 



provides the information of the luminosity distance and the redshift. The SN la samples lie 
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in the redshift region < z < 1.55. The luminosity distance-redshift relation is given by 
We use the CMB shift parameter measured by the five-year WMAP observation 



1090.04 



dz 



R = \ln^Hl I = 1.710 ± 0.019, (5) 







where Hq is the Hubble constant and is the dimensionless matter density at present. We 



use the BAO measurement from the joint analysis of the SDSS and 2dFGRS data 19|, |20| . 
which gives Dv{0.35)/ Dv{0.2) = 1.812 ± 0.060, where 

1/3 



-Dy(^BAo) - 

and Da{z) is the angular diameter distance 



(6) 



n / ^ 1 r dz' 



w-w' plan 



To obtain the constraints on the ly^'?// plane, we invoke a broadly used form of 
parametrization of the equation of state 



W{z) =Wq + Wail - a) =Wq + WaZj (1 + -Z) • (8) 

The constraint of ifo, Wa and fi^ is obtained by fitting the three parameters to this com- 
bined data set. The estimate of the parameters are found to be = — 0.891q;J4, Wa = 
— 0.18^q;74, fim = 0.25^o;q2- The two-dimensional constraint of WQ-Wa is obtained and 
shown in Fig. [2l 

B. Results of the constraints on the w—w' plane 

We reconstruct the the w-w' plane via Eq. ([8]) and 

W\z) = -aw a = -Wa/ {I + z) , (9) 

at late times for < 2; < 1. At each redshift, the two-dimensional constraint is obtained by 
converting the points on the boundaries of the confidence regions in the w^-Wa plane to the 
corresponding points in the w-w' plane, following Eq. ([8]) and Eq. Qj, for the 68.3% and 
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the 95.4% confidence regions respectively.^ Since the transformation between {w,w') and 
{wo,Wa) is linear, each point inside a confidence region in the WQ-Wa plane gives a distinct 
point inside the corresponding confidence region in the the w-w' plane. 

In the w-w' plane, we find that the cosmological constant, the phantom models, the 
up-rolling quintessence models, and the non-phantom barotropic fluids lie outside the 68.3% 
confidence region in the redshift regions < 2; < 1, 0.18 < z < 0.22, 0.4 < z < 1 and 0.7 < 
z < 1, respectively. This shows that the four models are ruled out at the 68.3% confidence 
level. On the contrary, the down-rolling quintessence models, including the thawing and the 
freezing models, overlap with the 68.3% confidence region for < 2; < 1. All of the models 
overlap with the 95.4% confidence region for < 2; < 1. Samples of the constraints on the 
w-w' plane at redshifts z = 0, 0.2, 0.7 and 1, together with the models, are shown in Fig. O 

IV. TEST OF THE METHOD 

In Sec. IIIIBt we invoked a criterion that a model is ruled out at the 68.3% confidence 
level if for some redshift the model's corresponding sector in the w-w' plane does not overlap 
with the confidence region at all. To test the validity of this criterion and to address the 
concern about the inherent bias of the parametrization against certain models, we perform 
a Monte Carlo test of our method. The criterion is invalid if for some redshift the resulting 
68.3% region from the Monte Carlo realization of the fiducial model does not overlap at all 
with the model's corresponding sector in the w-w' plane. We pick one or two fiducial models 
for each model category to test our method. 

The fiducial models used in the test include cosmological constant w{z) = —1, thawing: 
w{z) = -0.82 + 0.23 Ina + 0.08 (Zna)^ freezing: w{z) = -0.92 - 0.14 Ina - 0.05 {Ina^, 
up-rolling quintessence: w{z) = —1 + 0.0003 a~^'^, up-rolling phantom: (a) w{z) = —1.2, (b) 
w{z) = —1.16 — 0.2 Ina — 0.07 {InaY, down-rolling phantom: w{z) = —1 — 0.0003 a""^, and 
none-phantom barotropic fluids: w{z) = —1 + 0.0035 a~^. For all models, Qm is 0.25 and 
w{z > 1.55) is equal to w{z = 1.55). All models have w{z) < —0.8 for strong acceleration 

^ For the one-dimensional error propagation, following the reconstruction equations, the variance propaga- 
tion is Var{w) = {{w - (w))^) = {[wq - (wq) + (1 — a){wa - (wq))]^) = Var{wo) + (1 — a)'^Var{wa) + 
2(1 — a)Cov{wa,'Wa) and Var{w') = {—a)'^Var{wa)- 
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today. The trajectories of the models in the w-w' plane are shown in Fig. H] in the redshift 
region < 2; < 1. 

In the test, we realize each fiducial model by simulating the SN la, the BAO and the CMB 
data assuming current data quality. 1000 sets of simulated data are generated and fitted 
with the three parameters Wq, Wa and VLm- The values of Wq and Wa are converted to w and 
w' via Eq. ([8]) and Eq. ([9]) at each redshift. The 683 and the 954 of the 1000 Monte Carlo 
realizations, representing the 68.3% and the 95.4% region, are selected via their chi-square 
from the fiducial model. 

As a result, the criterion passes the test for all the models, that is, for each model the 
corresponding sector overlap with the 68.3% Monte Carlo realized region in the w-w' plane 
in the redshift region < z < 1. Samples of the test results at redshifts 2; = 0, 0.4, 0.7 
and 1 are shown in Fig. [51 It is seen that applying the criterion up to z = 1 might be 
pushing it to the limit, especially for the up-rolling quintessence, down-rolling phantom and 
the non-phantom barotropic fiuid cases. Yet the conclusion that the three models are ruled 
out at the 68.3% confidence level in Sec. IIII Bl still holds if we apply the the criterion only 
for < z < 0.7. 



V. CONCLUSION AND DISCUSSION 



Applying the bounds for various dark energy models in the w-w' plane for redshift < 
z < 1, we find that several models including the cosmological constant, phantom, non- 
phantom barotropic fiuids, and monotonic up-rolling quintessence are ruled out at the 68.3% 
confidence level based on the current observational data. On the other hand, down-rolling 
quintessence, including the thawing and the freezing models, is consistent with the current 
observations. All the models are still consistent with the data at the 95.4% confidence level. 



Using the same SN la data set, Shafieloo et al. |30| and Huang et al. |3l| also found the 
cosmological constant inconsistent with the data at the 68.3% confidence level. Barger et 
al. [2l| found the non-phantom barotropic fiuids excluded at the 95.4% confidence level based 



on the earlier data set. We notice that there was a time the observations favored w{z = Q) < 
— 1 jl3] but now the observations favor w{z = 0) > —1. However, the conclusions are drawn 
at the 68.3% confidence level at most. It is hoped that the next-generation observations will 
constrain the dark energy equation of state an order of magnitude better l|, ISj. We shall be 



7 



able to identify dark energy at higher confidence in the coming future. 
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< z < 1 
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W 



FIG. 1: Classification of dark energy models in the w-w' plane. Models are separated by the solid 
curves. The symbols "T", "F", "B", and "P" denote the "thawing", "freezing", "non-phantom 
barotropic", and "phantom" models, respectively. The quintessence models correspond to the 
region for w > —\ . The cosmological constant corresponds to the point (—1,0). The bold 
solid curve is both the lower bound for the non-phantom barotropic models and the bound that 
separates the down-rolling and up-rolling quintessence models (down-rolling takes the upper side) . 
The dotted curve is the bound that separates the down-rolling and up-rolling phantom models 
(up-rolling takes the lower side). 
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FIG. 2: The two-dimensional eonstraint of WQ-Wa based on the combined data set inchiding tlie 
Constitution set of SN la data, the CMB measurement from the five-year WMAP, and the BAO 
measurement from the SDSS and 2dFGRS. The dark and the hght gray areas correspond to the 
68.3% and the 95.4% confidence regions, respectively. The black point denotes the best-fit values. 
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FIG. 3: Samples of the constraints on the w-w' plane at redshifts z = 0, 0.2, 0.7 and 1. The dark 
and the light gray areas correspond to the 68.3% and the 95.4% confidence regions, respectively. See 
the caption in Fig. [T] for the description of the regions to that the models belong. The cosmological 
constant lie outside the 68.3% confidence region at all of the redshifts. The down-rolling phantom 
models lie outside the 68.3% confidence region at z = and z = 0.2. All the phantom models lie 
outside the 68.3% confidence region at z = 0.2. Both the up-rolling quintessence models and the 
non-phantom barotropic fluids lie outside the 68.3% confidence region at z = 0.7 and z = 1. The 
down-rolling quintessence models including the thawing and the freezing models overlap with the 
68.3% confidence region at the four redshifts. AS of the models overlap with the 95.4% confidence 

rpcrinn at thf fmir rpH shifts 



0<z< 1 



-1.2 -1.1 -1.0 -0.9 -0.8 




FIG. 4: Trajectories of the fiducial models for the test in the redshift region < z < 1. See the 
text in Sec. IIVI for the description of the models. The trajectories are shown by the black points 
and the black curves. 
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Cosnuill^cal ConslanI z = z = 0.4 z = 0.7 z = 1 




FIG. 5: Samples of the Monte Carlo test results of the method. From the first to the eighth row 

are the results corresponding to the fiducial models of cosmological constant, thawing, freezing, 

up-rolling quintessence, up-rolling phantom (a) and (b), down-rolling phantom and non-phantom 

barotropic fluids, respectively. The first column is the WQ-Wa plane and the rest four are the w-w' 

plane at redshifts z = 0, 0.4 0.7 and 1. The 68.3% Monte Carlo realized region is represented by 

the black points, while the 95.4% by the black and the dark grey points. The corresponding sectors 
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of the fiducial models are filled with the light grey color. See Sec. IIVI for discussion. 



